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We consider the problem of building an arbitrary N x N real orthogonal operator using 
a finite set, S, of elementary quantum optics gates operating on m < N modes - the 
problem of universality of S' on TV modes. In particular, we focus on the universality 
problem of an ra-mode beamsplitter. Using methods of control theory and some prop¬ 
erties of rotations in three dimensions, we prove that any nontrivial real 2-mode and 
‘almost’ any nontrivial real 3-mode beamsplitter is universal on m > 3 modes. 
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1 Introduction 

Around twenty years ago Reck etal. m considered the problem of building an arbitrary 
N x TV unitary operator using the set of elementary quantum optics gates. They gave a 
recursive algorithm which transforms an TV x N unitary matrix into an arrangement of 2-mode 
beamsplitters, phase shifters, and mirrors. If one wants to implement any TV x N unitary, 
however, then one needs to have access to all possible beamsplitters and phase shifters. This 
makes the results of [16; not particularly useful in a real experimental setting. Recently, 
Bouland and Aaronson p] considered the same problem albeit with a finite set of optical gates 
available. They showed that actually any single 2-mode beamsplitter that nontrivially acts 
on two modes densely generates unitary transformations 17(3) (or orthogonal transformations 
0(3), in the real case). Combining their result with the arguments from [16] they concluded 
that any nontrivial 2-mode beamsplitter densely generates U(m ) for in > 3 modes. For 
example, a real 2-mode beamsplitter is given by one rotation angle 8 £ [0, 27t] . Having 
such a beamsplitter we let it operate on pairs of three available modes. This way we get 
operators which effectively mix 3 modes. The resulting operators are dense in 0(3). As 
a direct consequence, for building orthogonal (unitary) transformation one does not need 
tunable beamsplitters - those whose 9 can be changed. In fact any 2-mode beamsplitter with 
a fixed 9 <£ {0, 7t/2, 7r, 37t/2} is universal for generation of quantum linear optics. They left the 
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problem of classifying optical gates that act on three or more modes open. This kind of gates 
can be easily built experimentally using coupled optical waveguides mm- We also note that 
there are some interesting developments in models of fermionic linear optics dang. In this 
paper we address the open problem given in [T] using methods that are orthogonal to those 
used in pQ. In a generic case, our approach can be used to obtain the desired classification 
for any number of modes which is impossible using representation theory arguments of [T]. 
We also give another proof of a 2-mode real beamsplitter universality which is entirely based 
on our method and does not make any use of the results of mm- 

Mathematically, a real m-mode beamsplitter is represented by an orthogonal m x m matrix. 
Throughout this paper we will follow convention of Nielsen and Chuang m and assume that 
this matrix has a determinant equal to one (some authors use a convention with — 1 CD- 
Under this assumption the set of m-mode beamsplitters forms the group SO(m). 

We say that a finite set of beamsplitters S C SO{m) is universal on m modes if and only 
if it generates SO(m), i.e. any m -mode beamsplitter can be approximated by a sequence of 
elements form S (or their inverses) with an arbitrary precision. This definition is analogous 
to the one for quantum gate universality. For example the famous set consisting of H and T 
gates is universal for one qubit as it generates any operation from SU( 2) with an arbitrary 
precision m- In this paper we consider the problem of m-mode universality when the set 
S is constructed (in some natural way which we explain in the subsequent sections) from a 
single 2- or 3-mode beamsplitter and show that: 

1. Any nontrivial 2-mode beamsplitter is universal on m > 3 modes. 

2. Almost any nontrivial 3-mode beamsplitter is universal on m > 3 modes. 

We also make several interesting statements for m-mode beamsplitters for an arbitrary m. 
They concern universality on k > m provided universality on m modes. The method we 
use to obtain the result is a combination of the fundamental theorem of control theory (see 
Theorem[l]) and some algebraic properties of the rotation group in three dimensions. Recently, 
Bouland and Aaronson [I] gave a proof of universality of a 2-mode beamsplitter (also for the 
case when SO( 2) is replaced by 517(2)). Their approach is based on representation theory 
and classification of subgroups of 517(3). As the authors of |1] point out, such classification is 
missing starting from SU (5) and therefore their approach has clear limitations. The method 
presented in this paper is complementary to p] and attacks the problem from a different 
direction. It divides the problem into two. For O m € SO(m) with the spectrum er(O m ) = 
{e, e - ^ 1 ,..., e i<?im / 2 , }, where each fa is an irrational multiple of 7r it boils down to 

proving that some particular elements generate the special orthogonal Lie algebra. This can 
be fully handled. For O m with at least one fa = an, where a € Q the subtle techniques to 
show that the product of two group elements that have finite order can have infinite order are 
required. We discuss this by considering an example with two rotations about the x and z axes 
by rational angle 9 whose cos(0) is algebraic of degree 2 in Section I3TT1 The techniques used 
in this section are based on cyclotomic polynomials and they were used in the similar context 
in m- They, however, do not generalise easily for an arbitrary rational angle. Therefore in 
the general case we use the recent results of Conway, Radin and Sadun concerning products 
of rotations mmmm- 
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Our method reveals the importance of mode permutations for m -mode beamsplitters, 
m > 3. The central role is played by the set S(O m ) = {P*O m P a : a g S m } where P CT are 
to x to matrices that permute modes of the considered beamsplitter. In particular we show 
that in a great number of cases universality of an m-mode beamsplitter O m on k > m modes 
reduces to showing that the set S(O m ) is universal on m-modes (rather than on to + 1 modes). 
Moreover, we show that already on 3-modes there is a beamsplitter that is not universal on 
3 and 4 modes. It corresponds to what we call the trivial action of permutation group, that 
is, to S(0 3 ) = {Os, Oa" 1 }- 

The paper is organized as follows. In Section[2]we discuss general aspects of the method we 
use in this paper. Next, in Section[3]we prove universality of a nontrivial 2-mode beamsplitter 
on to > 3 modes. The proof is divided into two parts. The first one is an elegant Lie-algebra 
calculation. The second is showing that the product of two finite order orthogonal rotations 
is a rotation by an angle which is an irrational multiple of 7 r. In the subsequent section we 
discuss some aspects of beamsplitters operating on higher number of modes, introduce S(O m ) 
and prove general results concerning S(O m )- Finally in Section [5] the 3-mode beamsplitters 
are discussed in details. 

2 The method 

In this section we sketch the method we will use throughout the paper. 

Let G be a connected Lie group and g its Lie algebra. We say that elements C 

G generate G, if and only if the set 

< gi ,...,g k >:= {g k \ ■ g k \ ■ ■ ■ g k ™ : a* g {1,..., k} , k t g Z, n g N}, 

is dense in G that is G = < g 1 ,,..., gk >■ Similarly we say that subgroups {.Hi,..., Hk} of 
G generate G iff the set 

>:= {g k l ■ g k a 2 2 ■ ■ ■ g k ” : <n € {1,..., k} , g ai € H ai , k r e Z, n € N}, 

is dense in G. Finally, let S = {Xi,... ,Xk} C g be a subset of Lie algebra g. We say that 
S generates g iff any element X of g can be expressed as a linear combination of X t 's and 
arbitrarily nested commutators of W’s: 


X — OLiXi + y oiij [Xi, Xj] -(- ^ otijk [Xi, [Xj, Xfc]] + ... 

3 i,3 i,j, k 

The following theorem [8] will be of the great importance in this paper. 

Theorem 1. Let G be a connected Lie group and g its Lie algebra. G is generated by one- 
parameter subgroups {e tx : t £ M}, X £ S where S is a finite subset of g if and only if S 
generates g as a Lie algebra. 

The problem which we are going to deal with is the following one: 

Problem 1. Let A = {ai,..., a^} C G be a finite subset of G. We want to show that A 
generates G, that is, the group generated by A is dense in G. 

We make use of Theorem [T] to solve Problem [l] To this end we note that one can always 
write 


at = e Xi , Xi g g. 
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We consider two cases: 

1. Assume that for all i G {1,..., k} we have that < ai > is dense in {e tXi : t G R}, i.e. 
< a,i> = {e tXi : t G R}. Under this assumption, by Theorem[l]we get that A generates 
G if and only if {Xi, ..., X/.} generates g. 

2. Assume that for some at the group < at > is not dense in { e tXi : t G R} but 

{A'i,...,Xfc} generate g. In this case we cannot directly apply Theorem [T] What 
we can do however is to replace ai by some element bi = e' ti that belong to < A > and 
is such that < > is dense in {e tYi : t G R} and {Xi,..., Y),..., X*,} generate g. If 

this kind of manipulation can be done for each “bad” ai G A then the problem is solved 
by means of Theorem [TJ 

In the following we will show that this is the case for beamsplitters. 


3 Universality of a real 2-mode beamsplitter 

A 2-mode beamsplitter is given by a matrix gg G £>0(2) of the form 


( cos(0) sin(0) \ 

^ — sin(0) cos^d) J ’ 


where 9 G [0, 27r[. Let us first look at the spectrum of gg. The characteristic equation reads: 


A 2 — 2Acos(0) + 1 = 0. 


And therefore spectrum is given by {e l£> , e l6 }. 


We want to show that three matrices: 


cos(0) sin(0) 0 

— sin(0) cos( 

0 0 


0 

\ 

/ cos(0) 

0 

sin(0) 

0 

, Ol,3 

0 

1 

0 

1 

J 

y — sin(0) 

0 

cos(0) 

1 

0 

0 \ 



0 

cos(0) 

sin(0) 



0 

— sin(0) 

cos(0) J 




(1) 


for a given and fixed 9 G [0, 27r[\{0, ir, |7r} generate 5'0(3)|? Following the reasoning 
plained in Section [2] we note that Oi t j = e Xi ’ j , where 


ex- 


0 9 0 

Xl,2 =1-000 


OO0\ / 0 0 0 

Ad,3 =1 0 0 0, X 2 ,3 = 0 0 0 

-0 0 0 / \ o -0 o 


0 0 0 

It is obvious that {X 1>2 , Xi j3 ,X 2j 3 } generate Lie algebra so(3) iff 0 ^ 0, as matrices 



> -Ei,3 — 



l The excluded angles correspond to either permutation of modes or =b identity operation. 
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form a standard basis of so(3). We also note that for 9 = an where a is irrational we 
have that < Oij > is dense in {e tXij : t G R}. Therefore using Theorem Q] we have that 
{0i,2, 01,3, 02,3} generate S'0(3) for 6 = an with irrational a. We still need to examine the 
case when a £ Q that is the case when for all three groups < Oij > ^ {e tXiJ : t £ R}. To 
this end we choose three new matrices that belong to < 0i,2, Oi,3, 02,3 >, he.: 

01 , 201 , 3 ; 01,202 3 ; 01 302 3 ' 


Our goal is to show: 

Statement 1. The Lie algebra elements corresponding to {Oi,20i,3, 01,202,3, 01 , 302 , 3 } 
form a basis of so(3). 

Statement 2. The elements {Oi,20i,3, 01,202,3, 01 , 302 , 3 } are rotations by angles that are 
not rational multiples of n , or equivalently there is no power k £ N for which they 
become identity matrix. 

The proofs of these two statements are given in the next three subsections. 

3.1 Finite and infinite order elements - examples. 

In this section we consider examples showing that the product of two orthogonal rotations 
by an angle 9 which is a rational multiple of n can be a rotation by an angle a which is not 
a rational multiple of n. In particular we show that this is the case for all 0’s whose cos(0) 
is algebraic of degree two. Our approach makes a heavy use of cyclotomic polynomials as 
they are useful in showing that for a complex number e la its argument a in not a rational 
multiple of n. These techniques were also used in to show that gates H and T generate 
SU(2), where the product of two specific rotation was considered. Here we provide much more 
general discussion which reveals the natural limitations of this method. The main purpose of 
this section is to provide some hands on examples and explicit calculations. For an arbitrary 
9 we will use a different approach (see section IQ) . 

3.1.1 Cyclotomic ploynomials 

The method we use is based on some properties of cyclotomic polynomials and therefore we 
start with their discussion (see [5] for more details). 

Definition 1. The n-th cyclotomic polynomial for n £N is given by 

®n(x) = n { x - e2i,r ”) > 

1 < k < n 
gcd (fc, n) = 1 

where gcd(fc,n) is the greatest common divisor of n and k. 

Cyclotomic polynomials have several useful properties. In the following we will use three 
of them: (1) cyclotomic polynomials are nronic and irreducible over Q, (2) Coefficients of 
cyclotomic polynomials are integers, (3) For any q £ N we have 

x q - 1 = JJ$ d (a:)- 

d\q 


( 2 ) 
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The first property is actually a nontrivial result due to Gauss (it can be however easily proved 
for prime n using Eisenstein criterion for irreducibility) [5]. The second one is a direct result 
of 

x q -l = H (z- e 2 iw t), 

1 <k<q 

and the definition of a cyclotomic polynomial. For a G C we say that the monic irreducible 
polynomial m a (x ) G Q[x] is the minimal polynomial for a over Q iff m a (a) =0. If a G C 
has a minimal polynomial over Q then we call it algebraic. Otherwise it is transcendental. 
The algebraic degree of a is a degree of its minimal polynomial. Now we can state the main 
theorem (cf. lemma 3.4 of H2) : 

Theorem 2. The following two are equivalent: (1) a G Q, (2) the minimal polynomial for 
a = e l2jra over Q exists and is cyclotomic. 

Proof. If a = p/q then (e l27ra )« = 1 and therefore the minimal polynomial m a (x) exists as a 
satisfies x q — 1 = 0. By 0 we know that a is a root of some <&d{x) where d\q. But $d(x) is 
irreducible and monic hence it is the minimal polynomial for a over Q. Conversely, assume 
m a exists and is cyclotomic. Then we have 

0 = m a (a) = $„(a) = II d |„$ d (a) = a n - 1. (3) 

By 0 we get e 2 ™ cn = 1 and hence c G Q. □ 

3.1.2 Products of rotations 

Theorem [2] can be used in particular for showing that for a given complex number e za the 
angle a is not a rational multiple of 7r. In this case it is enough to prove that either the 
minimal polynomial does not exist or it exists and is not cyclotomic. In the following we use 
it to study the composition of two rotations about orthogonal axes by an angle 9 which is a 
rational multiple of i r. Let 0\,2 and 02,3 be as in 0. We have 

trOi j2 02,3 = 2 cos (6) + cos 2 (0). 

From the other hand 0i, 2 02,3 is a rotation by an angle a and hence tr0i i 20 2 ,3 = 1 + 2 cos(a). 
As a result we get the equation which relates 9 and a: 

2 cos(a) = 2 cos (9) + cos 2 (0) — 1. (4) 

We first determine if e la is algebraic or transcendental. 

Fact 1. For a given by 0 e la is an algebraic number. 

Proof. As 9 is a rational multiple of n, using De Moivre’s formula we get that cos (9) is an 
algebraic number. Next, it is known that the sum and product of two algebraic numbers 
is again algebraic number. Applying this to 0 we get that cos(a) is an algebraic number. 
A complex number is algebraic iff its both real and imaginary parts are algebraic. To show 
that sin(a) is algebraic note that the field extensions Q [cos(a)] : Q and Q [cos(a), sin(a)] : 
Q [cos(a)] are both algebraic and consequently by the chain rule Q [sin(a)] : Q is algebraic. 
Hence e la = cos a + i sin a is algebraic. □ 
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Note that by Fact l.'l. 1 .~2l e la is never transcendental so in order to use Theorem [2] we need to 
determine its minimal polynomial. Putting x = e la = cos(a) + iyj 1 — cos 2 (a) we get 

x 2 — 2 cos(a)x + 1 = 0, (5) 

where cos(a) is determined by cos(0). 

Fact 2. The algebraic degree of cos(a) divides the algebraic degree of cos (0) and the algebraic 
degree of e la is twice the algebraic degree of cos(a). 

Proof. It is easy to see that the field extensions satisfy Q[cos(a)] C Q[cos(0)]. Using the chain 
rule for fields HU we get 

[Q[cos(0)] : <Q>] = [Q[cos(0)] : Q[cos(a)][Q[cos(a)] : Q] (6) 

As [Q[cos(0)] : Q] and [Q[cos(a)] : Q] are algebraic degrees of cos(0) and cos(a) respectively 
we get the conclusion. The relation between algebraic degrees e la and cos(a) comes from 
equation ©. □ 

Note that cos(0) can have arbitrary large algebraic degree and therefore minimal poly¬ 

nomial for e lCL can have any order. In the following we consider example when degree of 
cos(0) is 2 which by Fact [2] means the algebraic degree of e*“ can be either 2 or 4. As by our 
assumption cos(0) = a + by/C we get cos(a) = A + ByfC , where A and B are given in terms 
of o, b , c. The minimal polynomial of e la is: 

x 4 - 4Ax 3 + (4 A 2 + 2)x 2 - 4 Ax - 4B 2 C + 1 = 0 (7) 

We next determine possible values of a and b. Using De Moivre’s formula one can easily see 

that roots aq and aq with |aq| < 2 and |x 2 | < 2 of 

x 2 + ax + 6 = 0, (8) 

where a,b £Z are only possible values of 2 cos(0) of algebraic degree two. Direct calculation 
leads to: (1) cos(0) = ±-^=, (2) cos(0) = ±^, (3) cos(0) = ±j ± The corresponding 
angles are: cos(7r/4) = l/\/2, cos(7r/3) = y/Z/2, cos(7r/5) = 1/4 + l/4-\/5, cos(27t/5) = 
-1/4+ 1/4A/5, cos(3tt/5) = 1/4- 1/4^5, cos(4tt/5) = -1/4- 1/4^, cos(3tt/4) = -l/>/2, 
cos(27t/3) = —\/3/2. One explicitly checks that in all case polynomial © is not cyclotomic. 
This way we showed that the product of two rotations about orthogonal axes by rational 
angle 0, whose cos 9 is algebraic of degree two, is a rotation by an angle which is an irrational 
multiple of n. 

3.2 The proof of Statement 1 

We will make use of a compact form of Baker-Campbell-Hausdorff (BCH) formula for the 
group SO( 3). For the detailed derivation see [B] (cf. |7]). Let us first recall the definition of 
the BCH formula for an arbitrary compact semisimple matrix Lie algebra g. For X,Y £ g we 
define BCH(X, Y) in the following way: 

e BCH(x,y) = e .Y e v 
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It is known that BCH(A', Y) is given by an infinite sum. In the case of so( 3 ), however, there 
is a particularly nice compact formula for BCH(X, Y). This is due to the following two facts: 


1 . For X € so( 3 ) the characteristic polynomial is given by p( A) = —A 3 — ||X|| 2 A, where 
||X|| 2 = ^tr(X 4 X). Therefore by the Cayley-Hamilton theorem X 3 = — 0 2 X, where 
0 = 11 * 11 - 


e* = I + X + ^X 2 + fx 3 


3!' 


... + — X r 


=I+X[ 1- 


02 

3! 


6> 4 9 6 

"BT "7T 

sin 9 


+ X Z 


n\ 

1 

2 ! ~ 


9 1 

4 ! 


= I 


9 


-X + 2 


sin 2 (0/2) 2 


0 2 


X . 


04 

6 ! 


( 9 ) 


2. Let O = e x be a rotation matrix from SO( 3) and let Z 
see that Z = t ^^-X. Therefore: 


0 2 ° . Using formula m we 


log(O) := X 


ll^ll 


Z. 


( 10 ) 


Having these the BCH formula for so( 3 ) can be easily calculated. One simply writes down 
the expression for e x using © and then calculates the logarithm using m- Details can 
be found in [6]. In the case where X and Y are orthogonal tr(X*Y) = 0 the formula reads: 

BCH(X, Y) = qX + f}Y + y[X, Y], (11) 


where 


sin 1 (d)a sin L (d) b sin x (d) c 

a =- 7 —-77, P = -7-T> 7 =-7—-777 


( 12 ) 


a = sin 0 cos 2 (</>/2), b = sin/> • cos 2 (0/2), c = - sin 0 sin <f>, 


d= y/a 2 + b 2 +c 2 , 0 = ||X||, <)>= ||Y||. 


We need to show that 


{BCH(X 1 j 2 , X li3 ), BCH(X 1 i2 ,X 2)3 ), BCH(X 1)3 ,X 2)3 )} , (13) 

form a basis of so(3). Using (fill) and (TH21) we easily find: 

BCH(Xi i2 ,Xi j3 ) = S11 ^ 0 ^ sin s (^ cos2 ^* 1,2 + cos 2 ^Xi |3 - i sin 0X 2)3 ^ , (14) 

BCH(Xi i2 ,X 2j3 ) = S11 ^ ^ sin 0 ^cos 2 ^* 1,2 + cos 2 ~X 2 , 3 + ^ sin0X- li3 ^ , 

BCH(Xi i3 ,X 2)3 ) = S11 ^ ^ sin 0 ^cos 2 ^X 1>3 + cos 2 |x 2>3 - i sin0X lj2 ^ . 
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where d = sin#y 2cos 4 (#/2) + 1/4sin 2 9. Next we write down matrices (fill) in the standard 
basis of so (3) and get the following change of basis matrix 


sin 1 (d) . „ , sin 1 (d) . „ 

M = d Q sing - A so{3) = --—-sin 6 

The determinant of m is given by: 


(15) 


sin 1 (d) . 


sin# 


— 2 cos 6 (#/ 2 ) + - cos 4 (#/ 2 ) sin(#) + -sin 3 (#) 
2 8 


(16) 


To find its zeros of (fl6l) one can for example write down all functions in terms of t = tan(|) 
and solve the polynomial equations with respect to t.The relevant part of (1161) reads: 


— 2 cos 6 (#/ 2 ) + - cos 4 (#/ 2 ) sin(#) + -sin 3 (#) = 0 . 

2 8 

using sin 6 = 2 sin | cos | one gets: 

cos 3 (#/ 2 ) (— 2cos 3 (#/2) + cos 2 (#/ 2 ) sin(#/ 2 ) + sin 3 (#/ 2 )) = 0 . 


(17) 


(18) 


The first factor of (1T51) gives 9 = n. Using cos#/ 2 = and sin#/ 2 = A/j?, for the second 
factor of (1151) one gets polynomial equation that has only one positive real root t = y/2 — 1 (the 
remaining roots are t = —\[2 — 1 and four complex roots). That means # = 7r/2. Therefore 
the determinant m vanishes iff # = 0, # = n or # = tt/2. Hence (fl3l) form a basis of so(3) 
in all cases we are interested in. 

Remark 1. Using isomorphism SO( 3) = SU(2)/T.2 one can also prove this result working 
with SU{2 ) matrices. 


3.3 The proof of Statement 2 

For the proof we use the following result of m- 

Lemma 1 . Let A and B he rotations about orthogonal axes by 2-7r/p and 27r/g respectively. 
Consider the word: 

A ai B bl .■■A a ”-B bn , 

where none of ai’s are multiple of p/2 and none ofbi’s are multiple of q/2. If no two con¬ 
secutive terms represent rotations by n/2 and if the word is nonempty, then the word is not 
equal to identity. 

Corollary 1 . Let 0\,02 € SO( 3) be rotations about orthogonal axes by 2i r| and 27r^ ; re¬ 
spectively, where a < p and b < q and fractions a/p and b/q are not equal {1, |}. Then 

for any n£Na/e have (O1O2)" 7^ I ■ Thus the product of two finite order orthogonal rotations 
is a rotation by an angle which is not a rational multiple ofn. 

Proof. Let A and B be orthogonal rotations by ^ and =£■ respectively. We put 0\ = A a 
and O 2 = B b . Then (O 1 O 2 )" = ( A a B b ) n . By the assumptions a 7 ^ f and b 7 ^ |. Moreover, 
neither 0\ nor O 2 is a rotation by n/2. Therefore, using Lemma|I]we get ( A a B bs j n ^ 1. □ 
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Combining corollary [T| with the results of Section ITTT21 we get: 

Theorem 3. Assume 6 ^ {0, n, §7r}. Let 

S = {O k ,i ■ Ok,i £ 50(3), fc, l e { 1 , 2 , 3} , k < 1}, 

where Ok,i are given by W- The set S generates 50(3) and hence a real 2-mode beamsplitter 
is universal on 3-modes. 


3-4 A real 2-mode beamsplitter is universal on m modes, m > 3 

In the previous section we discussed the first non-trivial case, that is, 50(3). It turns out that 
having the result for 50(3) is almost enough to state the corresponding one for SO(N). To 
this end let us denote by {|f)}()Li a basis of C^. In this section we will consider the following 
set 


5 = {O m : Ok,i € SO(N), k,l G {1,.. - ,N} , k < 1} , (19) 

where Ok,i represent the matrix of a beamsplitter acting on modes fc, l: 

Oki = cos# (\k)(k\ + |^)(/|) + sin# (|fc)(/| - |Z)(fc|). 

The number of these matrices is N(N — l)/2. The Lie algebra elements satisfying e Xkl = Oki 
are given by 

X kl =6(\k)(l\-\l)(k\) = eE kl , (20) 

where {Eki}k<i is a standard basis of so(N) and we have the following commutation relations: 

[Xk,l , Xk,m] = #A^ m , [Xk,l > ^C,m] = b)Ak.mi = 

where we assumed k < l, l < to, k < m. 

We first, in section 13.4.11 give an argument based purely on Theorem Q] which shows that 
generation of 50(3) by {O 02 , 3 , Oi^} implies generation of SO(N) by 5 (defined as 
in (fT9l) ). This reasoning is then extended in section [4] to state some general facts about 
universality of m-mode beamsplitters, where to > 2. In section [3.4.2l we give another argument 
which is tailored specifically for 2-mode beamsplitters. It has an interesting advantage over 
the general argument which we discuss in section 13.4.11 

3.4-1 General argument 

Fact 3. Assume that { 0 ^ 2 , 02 , 3 , 0 1 , 3 } generates 50(3). Then 

{O k ,i : O k ,i e SO{N), k,l € {1,..., N} , k < 1} (21) 


generates SO(N). 

Proof. Assume that {Oi, 2 , 02 , 3 , 01 , 3 } generates 50(3). In particular, it means that after 
closure one can obtain any 2-mode beamsplitter O S 50(2) C 50(3) acting on any pair of 
available 3 modes. Therefore we have (at least in the limit) access to elements e Bl2 , e El3 , e E23 . 
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Repeating that argument for any three out of N modes we obtain all possible elements e EkI 
with 1 < k < l < N- some are obtained more than once. Note that set {E k {\ k< i is a standard 
basis of so(N) and the rotation angle of e Ekl is 1 which is clearly not a rational multiple of 
7 t. Therefore, by Theorem [lj we get the desired result. □ 


In the proof of Fact [3] we used that certain elements are available after closure. It is not 
clear, and in general may not be true, that they are available before closing set < S >. In 
the situation when the dense set generated by {Oi, 2 i 02 , 3 , Oi, 3 } does not contain elements 
e El2 , e Bl3 , e E23 the argument described above uses elements that are available only in the 
approximate sense to show generation of SO(N). From the mathematical point of view this 
is not a problem. One can also say that perhaps other elements e^ 12El2 , e^ 13El3 , e^ 23E23 with 
( f>ij not rational multiples of n are in fact available in < Oy 2 , 02,3, 0i,3 >• This can be true, 
however, it is not clear which ifiij's are possible and which are not. Therefore, the proof does 
not give any insight into how elements of SO(N ) are generated. Form the practical point 
of view one would like to know at least one example of a small number of elements that 
belong to < S >, that can be constructed in a simple way from the available beamsplitter 
and that enable generation of any element of SO(N). In section l3.4.2l we show an exemplary 
construction which provides the set of elements S^' 1 that generates a dense set in SO(N) and 
satisfies the assumptions of theorem [T] Moreover these elements are available in < S > before 
any closure and they are given by products of the original beamsplitter acting on selected 
triplets of modes. The construction can be viewed as an alternative proof of Fact [3] 

3-4-2 The example 

If 9 ^ 0 then {X k {\ k<l ,where X k i is given by CU1) . spans Lie algebra so(3). Moreover when 
6 = an, a (j Q then < O k ,i > is dense in { e tXk ■' : t G i?}. Therefore we can use Theorem 
|T]and obtain that {O k i} k<l generates dense subset of SO(N). What is left is to consider 
the case when a G Q. Note that for k < l and m < n , we have Ofe,;O m ,n = Om.nOfc.z iff 
{k,l} n {m,n} = 0. In this case non-trivial elements of the spectrum (i.e. those different 
from 1) of Ojt,iOm,n are nontrivial elements of spectra of O k ,i and O m>rl and hence if a is 
rational they are rational multiples of 7r. Thus we are interested only in the case when 
{k, 1} fl {to, n} ^ 0. Without any loss of generality we can assume that k < to. We have three 
possibilities: 


k = m : Ofc,/O fc ,n, 

l = m : O k jOi, n , (22) 

l — n O k jOm,i- 

Consider now the isomorphism: | k) ^ 11), |Z) >->• |2), |to) i-» |3)( | n) >->• |3)) between the 
3-dimensional spaces: Span c {|/c), |^), |TO)(|n))} and C 3 . Under this isomorphism we can apply 
Theorem Q] and obtain that spectra of matrices (l22l) are irrational multiples of n. Moreover, 
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we can use formulas m to find the corresponding BCH elements belonging to so (TV): 

BCH (X k ,i,X k>n ) = S11 ^. ^ sin e (^ cos2 +cos 2 ^ X k , n - isin OX^rfj , 

BCH (X k> i,Xi in ) = S11 ^ ^ sin 9 ^cos 2 ^X k>l + cos 2 ^X l>n + i sinOX^^ , 

HCH(A -k,rriT Xi ^ rn ') — ^ ^ sin 0 ^COS — V k,m T cos ~^X krn — sin , 

where d = sin d^j 2 cos 4 (0/2) + 1/4 sin 2 6. Next we explain how to chose indices {k,l} and 
{to, n} so that the corresponding BCH elements generate so (TV). This is in fact the nontrivial 
part of the extension from so (3) to so (AT), N > 3. 



Fig. 1. The pictorial representation of the rules for choosing elements Ofc,zO m , n . An index ( k,l ) 
represent the matrix i- If indices ( k , V) and (m, n) are connected by a line then matrix O & iOm,n 
is chosen as an element of a new generating set for SO(N). 


In case of SO( 3) we chosen as a new generating set 

S (3) = {O12O13, O12O23; O13O23}; 

and showed that corresponding BCH elements form the basis of so (3). The construction of 
the basis for TV > 3 is presented pictorially in figure [Q Each “box” with an index (i,j) 
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represents the matrix Oij. If there is a line between two boxes (fc, l) and (m, n), where k < m 
the product Ofc,zO m , n is a member of a new generating set for SO(N). For example in case 
of 50(4) we have: 


5( 4 ) 


{02,302,4, 02,303,4, 02,403,4} U {Oi,202,3, Oi,202,4, Oi, 302 , 3 } — 
~ 5 (3) UP (4) , 


where = {Oi,20 2 ,3, Oi, 2 0 2 > 4 , Oi, 302 , 3 }- The first set, 5 (3 \ is isomorphic to the one we 
used for 50(3) (under isomorphism |1) |2), |2) i-» |3), |3) i-G |4)). We need to show that 

BCH elements 


{BCH(X 2 , 3 * 2 , 4 ), BCH(X 2 , 3 X 3 , 4 ), BCH(X 2 , 4 X 3 , 4 )} U 
U{BCH(Xi, 2 X 2 , 3 ), BCH(Oi, 2 0 2 , 4 ), BCH(Oi, 3 0 2 , 3 )}, 

are linearly independent. To this end we write them in the standard basis {Ek,i} k<[ of so(4). 
The corresponding coefficients form columns of the change of basis matrix. We use the or¬ 
dering of the basis elements that reflects the structure of 5^ that is {^ 2 , 3 , p 2 , 4 , £ 3 , 4 } U 
{E12, £13, if 14}. Under these assumption the change of basis matrix has the following struc¬ 
ture: 


sin 1 (d) . . , sin 1 (d) . 

d .g sm O'Asow = 


sin 6» ( Aso ^ N 

d-e \ 0 3x3 Pso( 4) 


where 


( cos 2 f 

0 

cos 2 | \ 


( cos 2 f 

cos 2 f 

-2 sin 0 \ 

0 

cos 2 | 

0 

| , PsO( 4) — | 

\ sin# 

0 

COS 2 

V 0 

0 

0 ) 


^ 0 

| sin 6 

0 ) 


N = 


Since we are interested in the determinant only, the matrix N is irrelevant for calculation of 
det(A SO (4)) = det(A SO ( 3 )) • det(P SO (4))- But det(P SO (4)) = -§ sin# (\ sin 2 0 + cos 4 §) and 
hence vanishes iff 9 = 0 or 9 = n. 


General N To prove that BCH elements corresponding to products of rotations chosen 
according to the rules explained in figure 1 form the basis of so(N) we proceed by induction. 
Assume that det(A SO (jv-i)) is nontrivial except for 0 G {0, §, 7r, |7 t}. Next note that ~ 
5 (w-i) y r(n) w here 

R^) = {Ol,202,3, Ol,202,4, • • • , 0l,202,JV, Ol,302,3} • 


Therefore the change of basis matrix has a structure 


sin (d) . . sin 1 (d) . f A SO (n- 1 ) 


N 

PsO(N) 
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where 


PsO{N) 


( — \ sin 9 cos 2 | cos 2 f cos 2 § N 

cos 2 | ^ sin 6 0 0 

i sin 9 : 


V 


0 

\ sin 9 ) 


Therefore det Pso(N) = (— \) N sin^ 9 — 


cos 2 § det P' so(N) 


where 


( cos 2 | cos 2 | cos 2 | cos 2 | \ 

0 | sin 6 0 0 

P'so(n) = 0 0 |sinfi : 


V 0 0 


0 

\ sin 9 ) 


and is TV — 1 x N — 1 matrix. Hence 


det Pso(N) 



N—2 


■ N—2 n 4 0 

sin 9 cos - 
2 



N-2 


and det Pso(N) = 0iffd = 0ord = 7r which finishes the induction step of the proof. 
Theorem 4. Assume 9 ^ {0, n, A^-}. Let 

S = {O k>l : O k>l eSO(N), k,le{l,...,N},k<l}, 

where O k ,i represent the matrix of a beamsplitter acting on modes k,l: 

On = cos9 (\k)(k\ + |I)(Z|) + sind (|fc)(/| - \l)(k\). 

The set S generates SO(N). 

Remark 2. In \lfAj it was shown that one can not make arbitrary unitary transformations 
using only beam splitters when acting on two qubits. This result is not in contradiction with 
the fact that a 2-mode beamsplitter generates SO(4). The point is that in the setting of i'19f 
‘modes ’ are divided into two pairs and the beamsplitters can only act separately on these pairs. 
In our setting they are allowed to act on any pair of modes and therefore they can generate 
more than SO(2) x SO( 2) which happens to be entire SO(4:). 

4 m-mode beamsplitters 

In case of a 2-mode beamsplitter the freedom stemming from mode permutations was not 
significant since it was changing O 2 into Of 1 . For ?n-mode beamsplitters with m > 3 we have 
two situations which should be treated differently. 

Let Sym(m) be a permutation group of m elements, |Sym(m)| = to!. We consider an 
TO-mode beamsplitter O m = e Am £ SO(m), m > 3. Taking into account all possible mode 
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permutations the starting point is not a single beamsplitter but rather a set of beamsplitters 
given by matrices: 


S(O m ) ■= {PjO m Pa ■ o- £ Sym(m)}, 

where P a is an to x m permutation matrix corresponding to a £ Sym( to). Note that by 
Schur’s lemma S(O m ) has always at least two elements as the only permutation invariant 
matrix in SO{in) is the identity matrix. However, we still can have two cases: 


1. Trivial action of Sym(m): S(O m ) consists of O m and O rr f. 
permutation group is exactly as for a 2-mode beamsplitter. 


In this case the action of 


2. Nontrivial action of Sym(m): S(O m ) has at least two non-commuting elements. Com¬ 
bining this with a standard result of Kuranishi [9] saying that the set of pairs that 
generate a semisimple Lie group G is open and dense in G x G then is a good chance 
set S(O m ) generates SO(m). 

From now on we will say that O m is universal on k > m modes iff the set S(O m ) is universal 
on k > m modes. 

Remark 3. If a beamsplitter O m falls into the first case then it cannot be universal on m- 
modes if only m of them are at our disposal. On the other hand, if O m belongs to the second 
case it can be universal on m-modes without usage of any additional modes. 

Having X £ so (to) we can embed it into so (to + k) in ( m+k ) natural ways and similarly 
having O £ SO(m) we can embed it into SO(m + k ) in ( m +™) natural ways by choosing m 
out of to + k modes and letting X or O operate on them. For example in case when to = 3 
and k = 1 we have 


so(3) 5 o(4) : 


* * 

* * 

* * 


0 0 


* 

* 

* 


0 



0**\/00 0 0 \ 
000 || 0 ***| 

0 * * H 0 * * * 

0 * 0 / \ 0 * * * J 


50(3) ^ 50(4) : 


* * * 0 

* * * 0 

* * * 0 

0 0 0 1 


* * 

* * 

0 0 

* * 



We start with the following simple fact: 

Fact 4. Let X = {Xi,..., Xk} C so(to). Consider { m ^ k ) natural embedding of X into 
so(m + k). If X generates so(to ) then the { m ^ k ) natural embedding of X into so(to + k) 
generate so(m + k). 


Proof. As X generates so (to) we have a basis of so (m) at our disposal. In particular we can 
chose it to be a standard basis } );<? , where 1 < i, j < m. By definition of standard basis, 
it is obvious that natural embeddings of this basis into so(?n + k) gives a set which is 

an overcomplete basis of so (to + k). Therefore the result follows. □ 

Fact 5. The group < e Am > with A m £ so (to) is dense in { e tAm : t £ R} iff spectrum of 
e Am is given by : k £ {1,..., to}} with all f>i’s being irrational multiples of n. 


Proof. A direct generalization of 2-mode case. 


□ 
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Theorem 5. (Non-trivial action Sym(m) ) Assume that for an m-mode beamsplitter < S(Om) > 
is dense in SO(m). Then O m is universal on k modes for any k > m. 

Proof. As the set generated by < S(O m ) > is dense in SO{m ) we have access to any matrix 
from SO{m) (at least after closure). Having arbitrary SO(m ) matrix we choose the set 
B = {e Xi : Xi £ so(m), i £ {1,fc}}, where X = {Xi, ..., Xk} generate so(m) and spectra 
of Xfs are as in Fact[H By Fact 0] the ( m ^" fe ) natural embedding of X into so (m + k) generate 
so(m + k). As the assumptions of Fact [5] are satisfied we can use Theorem [T| and get the 
desired result. □ 

Theorem 6. (Trivial action ofSym(m)) Assume that for an m-mode beamsplitter S(O m ) = 

{ O m , 0~ j 1 } and the ( m (!" fe ) natural embeddings of S{O m ) into SO(m + k) generate dense set 
in SO{m + k) for some k. Then O m is universal on l modes for any l > m + k. 

Proof. The analogous to the proof of Theorem [5] □ 

Remark 4. (1) All the statements of this section remain true if we substitute SO(k), so(fc) 
with SU(k), su(fe). (2) Theorem [5| can be used to prove Theorem [|] We however found the 
calculation using only the BCH formula elegant and worth presenting. 

5 3-mode beamsplitter 

In this section we present the full discussion for S(0 3 ). 

The Lie Algebra so(3) is generated by {£ 12 , £ 13 , £ 23 }- Let O 3 = e As where 



aij £ R and ]£ \ctij\ 2 = 1. The permutation group consists of six elements 


Sym(3) = {e, (1,2), (1,3), (2,3), (1,2,3), (1,3,2)}. 


We want to consider set S(0 3 ) = {Pf0 3 P a : a G Sym(3)}. Making use of Pf0 3 P a = e p ° A3Pa , 
we have the following 5 (^ 3 ) set of Lie algebra elements 



Note that if S(0 3 ) = {0 3l 0 3 then PfA 3 P a = ±A 3 for any a £ Sym(3). It happens when 
ai 2 = —013 = a 2 3 , he. when 



In the following we divide our discussion of 3-mode beamsplitters into two cases according to 
the behavior under mode permutations. 
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5.1 Nontrivial action of Sym( 3 ) 

In this section we consider those beamsplitters whose A3 (E\ 2 — E13 + £23). Our goal 

is to first show that in such a case 5(^3) generate the Lie algebra so( 3 ). Let us start with 
the following simple lemma. 

Lemma 2 . Let X,Y € so( 3 ) be linearly independent. Then X, Y generate so( 3 ). 

Proof. Let {A12, Xl 3 , X23} be a basis of so( 3 ), where Xy = |*}(j| — |j)(i|. We have: 

X = OL\X x2 + P1X13 + 7 iX 23 , 

Y = 0:2X12 + P2X13 + 72X23 1 

[X, Y] = (71/?2 — / 3 i 72 ) X12 + (0772 — 710:2) X13 + {01&2 — 01^2) = 

= —M31X12 + M32X13 — M33X23 

We need to show that X , Y and [X, Y] form a basis of so( 3 ). To this end we calculate 

( 01 /?i 71 \ 

02 /?2 72 = - (M 3 2 i + M3 2 2 + M h ) ± 0 , 

—M31 M32 M33 ) 

since X and Y are linearly independent (at least one of My 7^ 0 ). 

Next we proceed with the proof that S(A 3 ) generates so( 3 ). There are a few cases to 
consider. 

1. Exactly one of aij is non-zero. Without loss of generality assume that 0,12 7^ 0 and 
oi 3 = 0 = 023. In this case elements 

A3 = GH2-E12, P£ 3 )A 3 P(i l3 ) = — 012 ^ 23 , -f > (2,3)^3-P(2,3) = ^ 12 ^ 13 , ( 24 ) 

form the basis of so( 3 ) and therefore matrices {Pj A 3 P a : a € S 3 } generate so( 3 ). 

2 . There are exactly two non-zero ay’s. Assume for example that a 3 2 = 0 . Then 

A 3 = 013^13 + 023^23: P(l } 3)-^-3-P(l,3) = — ^23^12 ~ ai3-El3j 

are clearly linearly independent and by Lemma [ 2 ] they generate so( 3 ). 

3 . All a t j \s are non-vanishing. We need to consider two cases: (a) a 3 2 — ai 3 7^ 023 

or ai2 = — a 13 7^ a 23 or a 12 = a 23 7^ — a 13 then A 3 and , 2 )'^ 3 -P(l 1 2) are linearly 
independent and by Lemma [ 5 ] they generate so( 3 ), (b) 0,12 7^ — ai 3 = a 2 3 then A3 and 
P(2 3) A3P^2.3) are linearly independent and by Lemma[ 2 ]they generate so( 3 ). 

Therefore we have 

Fact 6. Let O3 = e As be a beamsplitter which admits nontrivial action o/Sym( 3 ). Then 
5 (A3) generates so( 3 ). 


(23) 

□ 
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5.1.1 The case when 6 is an irrational multiple of n 
Combining Fact [6] with theorems [T| and [5] we obtain: 

Theorem 7. Let O 3 = e A3 be a 3-mode beamsplitter which admits nontrivial action o/Sym(3) 
and whose spectrum is given by {e lS , e~ lS , 1} where 9 is not a rational multiple of 7 r. Then 
5 ( 03 ) generates SO( 3) and O 3 is universal on k > 3 modes. 

Remark 5. Similar reasoning, however with more cases to consider, can be carried out for 
m-mode beamsplitters with m > 3. We will discuss generation of so (to) by S(A m ) for m > 3 
in a subsequent publication including analogous calculations for su(to) . 

5.1.2 The case when 9 is a rational multiple of n 

We start with a review of some important and relatively new facts concerning compositions 
of rotations in M 3 . 

Let 0\ and O 2 be two finite order rotations about axes separated by an angle a. In the 
series of papers mmmm Conway Radin and Sadun studied the group generated by 0\ 
and O 2 for large class of a’s. This group is characterized by relations that involve generators 
0\ and O 2 ■ In order to discuss these relations we recall some basic definitions from algebraic 
number theory. 

Definition 2. A complex number z G C is algebraic iff it is a root of some nonzero polynomial 
with rational coefficients. If z is not algebraic then it is called transcendental. 

Note that the set of algebraic numbers is countable. This can be easily inferred from the 
fact that there are countably many coefficients of polynomials in Q[a;] and each polynomial 
contributes finitely many algebraic numbers. Therefore the set of transcendental numbers is 
uncountable. 

Theorem 8. £1J/ Assume that for a G [0,27t[ there are nontrivial relations between 0\ and 
C? 2 - Then e l2a is algebraic. 

As an immediate consequence we get: 

Corollary 2. Assume e l2a is transcendental. Then there are no nontrivial relations between 
0\ and O 2 and in particular the group generated by 0\ and O 2 is infinite and dense in SO( 3) 
As e l2a is generically transcendental the above corollary covers almost all cases. What is 
left is the countable set of these a G [0, 27t[ for which e l2a is algebraic. Note that if a is a 
rational multiple of 7r then e 2la is algebraic as it is a root of unity (it satisfies x q — 1 = 0 for 
some q GN). 

Theorem 9. Let 0\ and O 2 be two finite order rotations about axes separated by an 
angle a which is a rational multiple of n. A group generated by 0\ and O 2 is infinite and 
dense in 50(3) with the following three exceptions: (a) Either 0\ = I or O 2 = I, (b) 0\ = I 
or O 2 = I and a = n/2, fcJOf = I = 0| . 

The exceptions to Theorem [9] correspond to Oi being rotations by: 0, n/2 , 7 r or 3n/2 
that we exclude. It is well known that when a is a rational multiple of 7r the algebraic order 
of e l2a can be arbitrary large. It is therefore interesting to understand better those angles 
for which e l2a has, for example, order two. Authors of [3] do this for the so-called geodetic 
angles, that is angles whose squared trigonometric functions are rational. For this kind of 
angle they prove relations between Oi and O 2 can occur only for finite number of a’s. The 
full understanding of all angles for which e l2a is algebraic is however still non-complete. 

Let us return to our problem. We assume O 3 = e Aa is nontrivial with respect to mode 
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permutations. Then the set S(Os) contains at least two rotations about axes separated by 
an angle a which is determined by coefficients of A3. Making use of the facts discussed above 
we have the following: 

Lemma 3 . Assume a is a rational multiple of it or is such that e l2a is transcendental. Then 
5 ( 6 ) 3 ) generates 50(3). 

Combining this with Fact [G] and Theorems Q] and [5] we get: 

Theorem 10 . Let O3 = e Aa be a beamsplitter which admits nontrivial action o/Sym(3) and 
whose spectrum is given by {e l6 , e~ l6 ,l} where 6 is a rational multiple of n. Let a be the 
angle between rotation axes of two different elements from 5(03). Assume a is a rational 
multiple of n or is such that e l2a is transcendental. Then O 3 is universal on k > 3 modes. 


5.2 Trivial action of Sym(3) 

In this section we show that when 5 ( 03 ) = { 03 , O^ 1 } that is for O 3 = e Az with 


A3 



12 ^ 


E 13 + E23 ), 


the group generated by the four natural embedding of O 3 into 50(4) is exactly 50(3). There¬ 
fore the beamsplitter given by O 3 is not universal on 3 or 4 modes. The four embedding of 
O 3 into 50(4) are given by O^fc = where 


4 ll 23 — E12 — E13 + E23, A234 — E23 — E24 + E34, 
A134 = E\3 — E\4 + E34, A124 = E\2 — E14 + E24. 


Elements Aijk are not linearly independent (A 123 + A 134 = A 124 + ^ 234 ) and one can easily 
verify that they span the 3-dimensional subspace 

Span R {Ti23, A234, ^4 134, 4I124} C so( 4 ). 

Lemma 4 . The space Span R {yli23, A234, A134, ^124} is a 3 -dimensional Lie subalgebra of 
so(4). 


Proof. It is enough to show that Span R {A 12 3 , T 234 , T 134 , A 12 4 } is closed under Lie bracket. 
The result follows from the commutations relations 


[A123, A234] = A134 + A124, [A123, ^134] = 4li24 ~ ^234, 
[A123, A124] = 4^234 ~ ^134, [^-234, ^134] = ^123 + A124, 

[A234, A124] = A123 — A134, [T134, A124] = A123 + A234. 


□ 

Let us remind that so(4) ~ so(3) © so(3). It is therefore natural to suspect that 

Span R {kli23, T234, T134, A124} ~ so(3). (25) 

Lemma 5 . The space Span R {kli23, ^234, A134, T124} is isomorphic to so(3). 
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Proof. Let 

X = — (^123 + A234 + A134 + A124) , 

y = ~ (dll23 + ^234 — ^134 — -A124) , 

Z = — — (A123 — ^234 — A134 + A124) • 

It is easy to verify that Span R {X, Y, Z} = Span R {A 12 3, T234, ^i34i ^124}- On the other 
hand we have: 


[X,Y]=Z, [Z,X]=Y, [V,Z\=X, 
which are commutation relations of so( 3 ). 


□ 


Combining Lemma [ 5 ] with Theorem Q] we get 

Theorem 11 . Let 9 be an irrational multiple of n. The group generated by four natural 
embeddings of O3 = e A3 into 50 ( 4 ) where 


A 3 



12 ^ 


E 13 + E23 ), 


is isomorphic to 50 ( 3 ) and therefore O3 is not universal on 4 modes. 

We are left with the case of 9 which is a rational multiple of n. By Lemma [ 5 ] the rotation 
matrices Oyfc act on some 3 -dimensional subspace of M 4 = Span R {ei, e2, e3, 64}. Note that, 
for example, O123 and O234 are rotations by 9 about axes given by: 


n 123 = 


n-234 = 


1 

7 ! 

1 


ei 


1 1 
e 2 - 7^ e 3, 


73 

1 


73 

1 


v ^' 62 ' 


e3 


= e 4 . 


Let a be the angle between 74123 and 74234 - One has cos(a) = — | and therefore a is an 
geodetic angle - an angle whose squared trigonometric functions are rational. The Primordial 
Theorem (Theorem 2 of [ 3 ]) lists all geodetic angles a = siiW 1 with p and q coprime 

that support nontrivial relations between two finite order rotations about axes separated by 
a. In our case sin(a) = ^/| and it does not belong to the list given in [ 3 ]. Therefore 
Theorem 12 . Let 9 be a rational multiple ofn. The group generated by four natural embed¬ 
dings of O3 = e As into 50 ( 4 ) where 


A3 



Ei 3 + E23 ), 


is isomorphic to 50 ( 3 ) and therefore O3 is not universal on 4 -modes. 

Remark 6. One can also show that O3 is not universal on 5 -modes - the natural ten em¬ 
beddings of O3 into 50 ( 5 ) generate the group isomorphic to 50 ( 4 ). We conjecture that on 
k-modes where k> 4 the (3) natural embeddings of O3 generate SO{k — 1 ). 
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6 Summary and outlook 

In this paper we discussed the universality problem of m-mode real beamsplitters for m = 2, 3 
from the perspective of control theory using some nice properties of the SO( 3) group. We 
also pointed out the importance of the set S(O m ) which is the orbit of adjoint action of 
permutation group through O m . In particular we showed that when S(0 3) = {O3, O^ 1 } 
the beamsplitter O 3 is not universal on both 3 and 4 modes (we also know it is the case 
for 5-modes). The study of similar phenomena in higher dimensions is a natural direction 
we want to explore. The other problem would be extension of the result presented here to 
complex beamsplitters. This requires proving several nontrivial results that we plan to discuss 
elsewhere. 
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